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Abstract 

We study forced oscillations of a viscoelastic rod with a body attached to its free end. 
The other end of the rod is fixed. Detailed analysis is presented for the case when the mass 
of the rod is negligible in comparison with the mass of the body. We assume a general 
form of distributed-order fractional constitutive equation for the rod. The existence of the 
solution for displacement and stress is proved. Some previous results are shown to be spacial 
cases of the present analysis. 

Keywords: fractional derivative, distributed-order fractional derivative, fractional vis- 
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1 Introduction 

In this paper we investigate forced oscillations of a body attached to a viscoelastic rod in the case 
when the mass of a rod is negligible in comparison to the mass of a body. A rod-body system 
is shown in Figure [TJ We consider the constitutive equation of a rod in a general form, given by 



u(x,f) 




Figure 1: System rod-body. 
QJjfy below, and investigate the function 
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obtained after applying the Laplace transform to the constitutive equation ([3]), connecting (j) a 
and 4> e which correspond to stress and strain. In this way we continue our approach studied in 
[5J U] in the context of solid-like rods. We find an explicit form of the solution as a convolution 
of a forcing function and a solution kernel for the case when the mass of a rod is negligible with 
respect to the mass of a body. 



1.1 Model 

First, we present equations in the case when the masses are comparable. Let a body of mass m 
be in translatory motion along the coordinate axis x that coincides with the rod axis. The length 
of the rod in undcformed state is L and its axis, at the initial time moment as well as during 
the motion, coincides with the x axis, see Figure [TJ Let x denote a position of a material point 
of the rod at the initial time to — 0. The position of this point at the time t > is x + u (x, t) . 
Then equations of motion of the rod-body system are 

d d 2 d 

—a(x,t) = p—u(x,t), e (x, t) = — u (x, t) , x e [0, L] , t > 0, (2) 
ox ot^ ox 

1 ,1 



4> a (n) D]a(x,t)d 1 = E / e (7) oD7eOM)d 7 , xe[0,L], t > 0, (3) 
o Jo 
d 

u(x,0)=0, —u(x,0)=0, <7(x,0)=0, e(x,0)=0, x&[0,L], (4) 

u(Q,t)=0, -Aa(L,t)+F(t) =m—u(L,t), t > 0. (5) 

ot z 

We use symbols a, u and e, in the equation of motion @i and in the strain ©2, to denote 
stress, displacement and strain, respectively, depending on the initial position x at time t, while 
p denotes the density of a material. Constitutive equation ^ corresponds to the distributed- 
order fractional derivative model of a viscoelastic body, E is a generalized Young modulus (a 
positive constant having the dimension of a stress) , <j) a and <p £ are given constitutive functions (or 
distributions), oDj is the left Riemann-Liouville fractional derivative operator of order 7 6 (0, 1) 
(see [22]) 

oVJy (t) := A (jr-^y * y (*)) , t > 0, (6) 

where T is the Euler gamma function and * is the convolution. Recall, if /, g G L] oc (R) , 
supp /, g C [0, 00) , then (/ * g) (t) := f* f (r) g (t - t) dr, tel. We refer to [H El |22] for 
the basic definitions and assertions of fractional calculus. Initial conditions in (jj) specify that 
the rod-body system is unstressed and in the state of rest at the initial time instant. Boundary 
condition ©i means that one end of the rod is fixed. The other boundary condition ([5J2 is the 
equation of translatory motion along the x axis of the body attached to the free end of the rod. 
In A stands for the cross-section area of the rod and F stands for the known external force 
acting on the body. 

Regarding the constitutive equation ([3]) we list some special cases. 

I Fractional Zener model of a viscoelastic body 

(l + a Q D?)<r(x,t) = E(l + b D?)e(x,t) (7) 

is obtained from (J3J by choosing 

^( 7 ):=5( 7 ) + a5(7-a), cj> e (7) := 8 (7) + bS (7 - a) , a G (0, 1) , < a < b, (8) 

where 5 denotes the Dirac delta distribution. 
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II Distributed-order model of a viscoelastic body 



f a 1 D]a(x,t)d-/ = E f V B]e (x,t) d 7 
Jo Jo 



(9) 



is obtained from ([3]) by choosing 

<$> a (7) := a\ 4> e (7) := b\ 7 6 (0, 1) , < a < 6. (10) 
We note that (|10[) is the simplest form of <p a and <p e providing dimensional homogeneity. 

We note that another area of application of fractional calculus in mechanics is related to non- 
local effect in elasticity. In this case derivatives of non-integer order are taken with respect to 
spacial coordinate. For this kind of applications see [151 116) . However, our constitutive equation 
© is of more general form than those used in [T3 [TB] ■ 

After giving a formal convolution of a forcing function and solution kernel, we analyze as- 
sumptions on M and obtain explicit solutions for the displacement u and stress a. This is done 
in the main Theorem [5] We refer to [2] for the analysis of oscillations of an elastic rod with the 
mass attached to its end. We also refer to [3J , where we treated a problem, similar to the present 
one, for a rod described by a fractional distributed-order model. 

We analyzed in [5J [S] the initial-boundary value problem © - © in the general framework, 
as well as for the cases of constitutive equation © given by ([7]) and ©. In [6l [7] the wave 
propagation in a viscoelastic solid-like rod (the constitutive equation was ©) of finite length 
with one of its ends fixed to a rigid wall. We considered two types of boundary conditions: 
a prescribed displacement and prescribed stress on rod's free end. Similar problem of wave 
propagation was analyzed in [3] for the fluid-like material, i.e., when the constitutive function is 
given by 

(1 + a D?) a (x, t) = E (b Df + h d£ 1 + b 2 Df 2 ) e (x, t) , 

where a, b , 61, b 2 are positive constants, < a < (3 < /3 1 < (3 2 < 1- We refer to [I31IIS1HI] for 
the detailed account of applications of fractional calculus in viscoelasticity. Problems similar to 
© - ([5]) were also treated in [THl HO] with the constitutive equations related to the distributed- 
order model © in the special cases. 

Note that our results are obtained by the use of the general framework including methods 
of the distribution theory. We refer to [33] for the definition of space of tempered distributions 
supported by [0, 00) , the Laplace transform within this space, as well as of the convolution. 

2 Convolution form of solutions 

The system © - © transforms into 

d d 2 d 

—a(x,t) = K 2 —u(x,t), e(x,t) = — u(x,t), xe[0,l],*>0, (11) 

1 ,1 

<M7)oDMM)d 7 = / &(7)oD?e(*>*)<17, *e[0,l], t>0, (12) 
Jo 

d 

u(x,0) = Q, —u(x,0) = 0, a(x,0) = 0, e(x,0) = 0, xg[0,1], (13) 
u(0,t) = 0, -<7(l,t) + F(t) = J^u(M), t>0. (14) 
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This is done by introducing the square root of the ratio between the masses of a rod and a body 



pAL 



and dimensionless quantities 

V AE vV AE ) vV AE ) 

In writing (fTTj) - ([T4"f we omitted bar over dimensionless quantities. Note that the choice of 
dimensionless quantities implies that the case of a rod without the attached mass (m = 0) 
cannot be studied as a special case of equations ([TT|) - (114)) . 

Now, we turn to the case of a light rod, i.e., the case when p = (n = 0). Then, the 
initial-boundary value problem (jlll) - (fT4"| takes the form 

d d 

— a(x,t)=0, e(x,t) = — u(x,t), if O,lU>0, (15) 
ox ox 



f ^( 7 ) D7a(^t)d 7 = / B ( 7 ) o D7e(a;,t)d7 ) xG [0,1], i > 0, 
o Jo 



(16) 

Jo 

u(a;,0) = 0, —u(x,0) = 0, a(x,0) = 0, e(x,0) = 0, xg[0,1], (17) 
u(0,t) = 0, -a (l,t) + F (i) = — u(l,t), t>0. (18) 

From (IT5|) i . we have that stress does not depend on the spatial coordinate, i.e., 

a = a (t) , t> 0. 

Then, ([TBI implies that the strain also does not depend on the spatial coordinate, since the 
left-hand side of (Tl6l) is not a function of x. Therefore, 



e = e (t) , t > 0, 

and from (1151) 9 we conclude that 

u(x,t) = xe{t) , a; G [0,1] , * > 0, (19) 

since u satisfies (I18lh . Taking into account these considerations, the initial-boundary value prob- 
lem (JT5J) - (3D becomes 

«^( 7 ) D7<7(*)d 7 = / e (7)oDMt)d7, *>0, (20) 
Jo 

- <r (t)+ J F(t) = ^j e (t), t>0, (21) 
«r(0) = 0, e(0)=0, ^(0)=0. (22) 

The initial-value problem (|2U|) - (|2"2")) has been studied in [3] in the case of the constitutive 
functions of the form (ITUl) . i.e., 

<j) a ( 7 ) = a 7 and </> £ ( 7 ) = cS ( 7 ) + 6 7 , < a < 6, c> 0. 
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We see that the motion of a system rod-body is determined only by the constitutive equation 
for the rod (j2"0)) and the equation of motion for the body ([2~Tj) . From (Ti"T)| it follows that e is the 
displacement of the body (and of a free end of a rod). 

In order to solve the system (|20|) . (j2Tj) subject to initial data (|22|l . we use the Laplace 
transform method. Recall, the Laplace transform of / £ L ; 1 oc (K) , / = in (—00, 0] and \f (t)\ < 
ce kt , t > 0, for some k > 0, is defined by 

/>oo 

/ (s) = C [/ (<)] (s) := / f(t) e~ st dt, Res>k 

JQ 

and analytically continued into the appropriate domain D. Recall, <S' (R) is the space of tempered 
distributions on R, while <S+ is its subspace consisting of tempered distributions supported by 
[0, 00) . We refer to [22] for the properties of this space as well as for the Laplace transform within 
it. We also use C ([0, 1] , <S>+) to denote the space of continuous functions on [0, 1] with the values 
in S' + . 

Applying formally the Laplace transform to (|2U| - (|2"21 we obtain 

u (s) [ 4> a (7) s 7 d 7 = e (s) [ £ (7) S ^d 7 , s S A (23) 



By we have 



a (s) + s 2 e (s) = F (s) , s E D. (24) 



&{x,s) = A ,p( s ^ ( x > s )i se -°' ( 25 ) 



where we introduced M by (TTJ). Note that for s 



E(w)=E (u) + iE (w) := = » , w G (0, 00) . (26) 

(iw) J Q 0^ (7) (iw) 7 d 7 

is a complex modulus, see [TU]. Functions E' and £"' are real- valued. They are called storage 
and loss modulus. 

By (HH an d (ESJ) we have 

e(«) = F(a)P(s), er(a)=J , (a)Q(s), s e D, (27) 

where 

P(s):= ^—5, Q(a):= sffl, (28) 

l + (sAf (,s)) 2 W l + (sM(s)) 2 

and M is given by (fT|). 

Inverting the Laplace transform in (|27[) . we obtain 

e(t) = F(t) *P(i), ir(t) = F(t)*Q(t), t>0. (29) 

We now discuss the restrictions that (f) a and <f> e must satisfy. There are two conditions that 
those functions must satisfy. First, P and Q must be real- valued functions, so that strain e 
and stress a, given by (|2"9"|) . are real. The second condition imposes the Second Law of Ther- 
modynamics, which requires that (in isothermal case) the dissipation work must be positive. 
Mathematically, these conditions read as follows. 

Condition 1 
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(i) There exists xq £ K such that 

Jo <$>a (7) ^7 



£ K, for all x > xq. 



M(x) 

(ii) For all u) 6 (0, 00) we have 

E' (lu) > 0, E" (oj) > 0, 
where E' and E" are storage and loss moduli, respectively, given by i26}) , see U Uf . 

The motivation for (i) of Condition Q] follows from the following theorem of Doetsch. 

Theorem 2 ([11, p. 293, Satz 2]) Let f(s) = C[F}(s), Res > x Q £ R, be real-valued on the 
real half-line s 6 (xo, 00) . Then function F is real-valued almost everywhere. 

Alternatively, if f is real-valued at a sequence of equidistant points on the real axis, then 
function F is real-valued almost everywhere. 

If 4> a and (f> £ are such that (i) of Condition [T] is satisfied, Theorem [2] ensures that inversions 
of (j27|) with ((28)) . given by (|29|) . are real. As it is well-known, [10], (ii) of Condition Q] guarantees 
that the Second Law of Thermodynamics for the isothermal process is satisfied. We shall see from 
Proposition[S] below that Condition [T] along with an additional assumption on the asymptotics of 
M (Assumption^] below) guarantees that the poles of the solution kernel in the Laplace domain 
([2"5)l belong to the left complex half-plane. In this case the amplitude of the solution decreases 
with the time; this is a characteristic behavior for a dissipative process. 

Remark 3 

(i) Condition Q] is satisfied for the fractional Zener model ^ and for the distributed- order 
model (0^. In the cases of ^ and (0), the function M has a form 



/ 1 _)_ as a 

M ^ = VTT6^' seC \(- TO '°]< 0<a<6, ae (0,1), (30) 



M(s) = y^^, S eC\(-oo ) 0] )0 <a< & , (31) 

respectively. Putting s — x € M in the previous expressions we see that (i) of Condition [7] 
is satisfied with xq — 0, while < a < b ensures that (ii) of Condition [7] is satisfied, see 

00123'. 

(ii) In general, stress a as a function of a real-valued strain s may not be real-valued. In J5}j 
we have such a situation, because (i) of Condition [JJ is not satisfied. This shows the role of 
Condition^ (i). 



3 Inversion of the Laplace transforms 

In order to obtain e and a by (|27[) , we have to obtain functions P and Q, i.e., to invert the 
Laplace transform in (1281) . Let us make an additional assumption on the function M, ([TJ. 



G 



Assumption 4 Let M be of the form 



M{s) 



r(s)+ih(s), as \s\ — > oo. 



Suppose that 



lim r (s) 



Coo > 0, lim ft, (s) = 0, lim M (s) = cq 



|s|— >oo s|— 



for some constants c^, cq > 0. 

Assumption 0] is inspired by the fractional Zener (JT]) and distributed-order model © . We 
note that both of these models describe the viscoelastic solid-like body. For the both models 



Then f has two different zeros: so and its complex conjugate so, located in the left complex 
half-plane (Res <0). The multiplicity of each zero is one. 

Proof. Since sM (s) = sM (s) (bar denotes the complex conjugation), it is clear from ([32]) that 
if so is zero of (|32p. then its complex conjugate sq is also. 

We use the argument principle in order to show that there are no zeros of (1321) in the upper 
right complex half-plane. Let us consider contour 7^ = j R1 U j R2 U 7^3 U 7^4, presented in 
Figure [5] Let j R1 : s = x, x £ [r, R] , where r is the radius of the inner quarter of circle and R 



is he radius of the outer quarter of circle. Then, by (i) of Condition [T] (with xq = 0) we have 
Re / (x) > and Im/ (x) = 0. Assumption 2] implies lim r _,.o f (x) = I and limfl_ ! . 00 / (x) = 00. 
Hence, Aarg/ (s) = for s 6 j m , as r — >• 0, R — >• 00. Let 7^3 : s = Re llp , ip £ [0, ^] . Then, by 
Assumption [H for ^ = we have 




(32) 



Tim s 




Figure 2: Integration contours j L and 7^. 



Re/(i?) 



c^R 2 -> 00 and Im/ (ii) -> as J? -» 00. 



For 



^ Assumption 2] implies 
Re/(ite f *) «-c : 



c^i? 2 -> -00 and Im / (ite 1 ^) -> as R ->• 00. 



(33) 
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By using Assumption 2] in (1321) we obtain 

Im / (Re iip ) « c^R 2 sin (2ip) > 0, <p G 



7T 

°'2 



Therefore, Aarg / (s) = tt for s G 7^2; as i? — >■ oo. Let 7^3 : s — cue 1 2 
have 

/ (iw) = 1 - lu 2 M 2 (uj) , cj e [r, i?] . 
Using 1(25)1 in ([33)) we obtain 



as i? -> 00. (34) 
iu>, uj G [r, J?] . Then we 
(35) 



/ (iw) = 1 - w 



(25») 2 + (i5"M) 2 



Thus, 



E" (uj) 

Im/(iw)=u; 2 5—^ , > 0, w G [r, i2] , 

(s'(w)) a + (£?"H) a 



(36) 



due to (it) of Condition [T] Moreover, Assumption 2] applied to ((33)) implies 

Re/(iw) w -cj^o; 2 ^ -00 and Im / (iu;) ->• 0, as R ->■ 00, (37) 
Re / (iw) 1 - c 2 w 2 -> 1 and Im / (ioj) ->■ 0, as r ->• 0. (38) 

We conclude that A arg / (s) = — 7r for s G 7fi3 , as r — ► 0, i? — >• 00. Let j m : s = re' v , tp G [0, J] . 
Assumption 21 for (p = and ^ = ? implies 

Re / (r) w 1 + c^r 2 -> 1 and Im / (rJ^Oasr-) 0, 

Re / (re 1 *) si - c 2 r 2 ->■ 1 and Im / (re 1 *) as r ->■ 0, (39) 

as well as 

Re / (re iv ) w 1 + cgr 2 cos (2<p) -> 1 and Im / (re iv ) « c^r 2 sin (293) -> 0, (40) 

for <p G [0, ^] , as r — > 0. We see that Aarg / (s) = for s G 7^4, as r — >• 0. Thus, we conclude 
that 

Aarg/(s) = for s G 7 fl , as r — > 0, i? — > 00. 

By the argument principle and that fact that the zeros are complex conjugated we conclude that 
/ has no zeros in the right complex half-plane. 

We shall again use the argument principle and show that there are two zeros of (j3"2"| in the 
left complex half-plane. Let us consider contour 7i = j L1 Uj L2 U7 i3 U7 i4 , presented in Figure 
BJ Let 7il : s = we 1 * = iw, w e [r, R] . Then p6 ]> . ([57 ] and ([38 ]) hold. Thus, A arg / (s) = tt for 
s G 7 L1 , as r -4 0, i? 00. Let 7 i2 : s = i?e' v , ^ G [f , ^f] . Then, for <p = | we have that ((33)) 
holds. For = 7r Assumption 2] implies 

Re / (i?e i7r ) « c^i? 2 ->■ 00 and Im / (i?e i7r ) as R -> 00, 



while for 9? = ^ it implies 
Re / (i?e i4 r) 
By d3U), we have 

Jmf(Re iv ) 

Im/fRe 1 ?) 



c^R -> —00 and Im / (iie 1 2 ) — > as i? — > 00 



c^sin(2^)<0 ) ^G 



37T 

c^i? 2 sin (2y?) >0, <^G (tt,— 



7r ] as R — > 00 . 



as R — >• 00. 
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We conclude that A arg / (s) = 2ir for s G 7 L2 , as i? — > cxd. Let 7 i3 : s = we 1 2 = — iu>, lu G [r, R] . 
By K32J) we have 



/ (-iw) = 1 - w'M 2 (iw), w G [r, i?] 



(41) 



since M (a) = M (s). Using (gU) in (gT]) we obtain 

2 £' (lu) + iE" (w 



Thus, 



/(-iw) = l-w' 



Im/(-iw) 



., , we [r, J?l . 

(B'( W )) 2 + (£"H) 2 



(E>(u>)) 2 + (E»(u)f 



< 0, w G [r,i?] 



due to (m) of Condition [T] Note that Assumption 2] applied to (|4"Tj) implies 



Re/(-iw) 
Re/(-iw) 



2 2 



— cxd and Im/ (iw) — >• 0, as R — > 00, 



1 - CqCj 2 — > 1 and Im / (icj) — >• 0, as r -» 0. 



We conclude that A arg / (s) = 7r for s G 7 i3 , as r — > 0, i? — s> cxd. Let j L4 : s = re lip , ¥>€[§, • 
Then, for ^3 = •? we have that (|39| holds. For y> = 4^ Assumption [4] implies 



Re 



/ (re 1 ^ j « 1 - Cpr 2 1 and Im / ^re' 32 " j -> as r -> 0. 



We have that also (f4*Uf holds for ^ 6 [|, 4^] , as r — » 0, so that A arg / (s) = for s G j Li , as 
r — > 0. Thus, the conclusion is that 

Aarg/(s) = 47r for s G j L , as r — » 0, R — > 00. 

This implies that / has two zeros in the left complex half-plane. ■ 

The following theorem is related to the existence of solutions to system (|20D - (|2"2"j) . 

Theorem 6 Let M satisfy Conditional^ with xq — and Assumption^ Let F G S' + . 

(i) The displacement u, given by H9\), as a part of the solution to \20jl - flUfy) , is given by 



and 
P(t) 



u (x, t) = xe (t) , where e (t) = F (t) * P (t) , x G [0, 1] , t > 
M 2 (qe- i7T ) 



- I Im 

7T 



1 + (qM (qe-^)Y 



e- 9t dq + 2Re(Res (P(s) 



e st ,s 



(42) 

t > 0, 
(43) 



P(t) = 0, t<0. 



The residue is given by 

Res (P(s)e st ,s a ) = 



M 2 (s) 



t > 0, 



(44) 



where f is given by L32\) and sq is the zero of f. Function P is real-valued, continuous on 
[0, cxd) and u G C ([0, 1] , <S+) • Moreover, if F is locally integrable on K. (and equals zero on 
( — 00, 0}) then u G C ([0, 1] x [0, 00)) . 



9 



(ii) The stress a as a part of the solution to H2U\) - is given by 

a(t) = F(t)*Q(t), xe[0,l],t> 0, 

where 



- I In, 



Q(t) 



Q(t) = 0, t<0. 
The residue is given by 



1 



1 + (qM (qe-™)y 



- 9t dq + 2Re(Res (q (s)e st , 



sq 



t > 0, 
(45) 



Res 



e st ,s 



t > 0, 



(46) 



where f is given by A Sty) and sq is the zero of f. Function Q is real-valued, continuous 
on [0,oo). Moreover, if F is locally integrable on R (and equals zero on (—00,0],) then 
a G C([0,oo)). 

Proof. First, we prove (i) . Condition Q] ensures that P is a real- valued function. We calculate 
P (t) , t e M, by the integration over the contour given in Figure |3l The Cauchy residues theorem 




\q S 



Figure 3: Integration contour T 



yields 



<j> P (s) e st ds = 2vri (Res (P (s) e st , s ) + Res (P (s) 



e s \s 



(47) 



where T = T x U T 2 U T 3 U T 4 U T 5 U T 6 U T 7 U -y , so that poles of P, given by (HHJ), lie inside the 
contour T. Proposition [S] implies that the pole sq, and its complex conjugate So, of P are simple. 
Then the residues in (|4"T)) can be calculated using (|44l) . 

Now, we calculate the integral over T in (|4"T|) . First, we consider the integral along contour 
ri = {s =p + iR\p e [0, s ] , R > 0} , where R is such that the poles so and so lie inside the 
contour T. By (|28l) and Assumption 2] we have 



P(s) 



< 



C 



(48) 
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Using (|48l) . we calculate the integral over Li as 



lim 

R^oo 



P (s) e st ds 



Vi 



< lim 

R— too 



,(p+ifl)t 



P + i/Z) 
< C lim / ° -^e p *dp = 0, i > 0. 



dp 



Similar arguments are valid for the integral along contour L7: 



lim 

R— >oo 



P (s) e s *ds 



0, t > 0. 



Next, we consider the integral along contour L2. By (l4"5f we obtain 



lim 

R— ¥00 



P (s) e st ds 



< lim 



< C lim 

-R— >oo 



P (Re 1 *) 
* 1 



i?£ cos 



P 



|iPe I0 |d0 
0, i > 0, 



since cos <j> < for <f> 6 [f ,?r] . Similarly, we have 



lim 

R— >oc 



P (s) e s *ds 



0, i > 0. 



Consider the integral along L4. Let \s\ — » 0. Then, by Assumption [H Af (s) — » Co and 
sM (s) -> 0. Hence, from ([251) we have 



P(s) « |M(s)| 2 w c§, as |s| -> 0. 
The integration along contour L4 gives 



(49) 



lim 

r->0 



P (s) e s *ds 



< lim 

r->0 



Plre 14 



lire 



r->0 



< eg lim / re r * cos <M<?(> = 0, t > 0, 



where we used ([41?]) . 

Integrals along L 3 , L 5 and 7 give (t > 0) 



lim / P(s)e st ds 

R— >oo / r 



lim / P(s)e st ds 



M 2 (qe^ 



1 + [qM (qe^)y 
M 2 (qe' iyT ) 



- qt dq, 



1 + {qM{qe-' m )) 



~ qt dq, 



lim / P(s)e st ds = 2vriP (t) . 



(50) 

(51) 
(52) 



We note that [[52]) is valid if the inversion of the Laplace transform exists, which is true since 
all the singularities of P are left from the line 7 and the estimates on P over j imply the 
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convergence of the integral. Summing up (1501) , (l5"Tj) and (1521 we obtain the left hand side of (jTT)) 
and finally P in the form given by (1431) . Analyzing separately 



Im 



M 2 (ge~ i7r ) 
1 + (qM (ge" i7r )) 2 



e _9t dg, 2 ]T Re (Res (P (s) e st , 



, s 



we conclude that both terms appearing in (l43l are continuous functions on t G [0, oo) . This 
implies that u is a continuous function on [0, 1] x [0, oo) . From the uniqueness of the Laplace 
transform it follows that u is unique. Since F belongs to <S^_, it follows that 

u(x,-) = x(F(-)*P(-))eS' + , 

for every x G [0, 1] and ueC ([0, 1] , S' + ) . Moreover, if F G Lj oc ([0, oo)) , then u G C ([0, 1] x [0, oo)) 
since P is continuous. 

Second, we prove (ii) . Again, Condition [T] ensures that Q is a real-valued function. We 
calculate Q (t) , t G R, by the integration over the same contour from Figure [3l The Cauchy 
residues theorem yields 



j> Q (s) e st ds = 2vri (Res (q (s) e st , s ) + Res (q (s) e st , 



■so 



(53) 



so that poles of Q lie inside the contour T. The poles sq and Sq of Q, given by (|2"8]l are the same 
as for the function P. Since the poles sq and sq are simple, the residues in (|53[) can be calculated 
using (|46l) . 

Let us calculate the integral over T in (|53|) . Consider the integral along contour 

r x = {s = p + m | p g [o, s ] , r > o} . 

By (|28p and Assumption [H we have 



Q(«) < 



C 



(54) 



Using (|54l) we calculate the integral over Ti as 



lim 

R— ¥00 



Q (s) e s *ds 



< lim 

-R— too 



,{p+iR)t 



Q(p + ii?) 

rso i 

< C lim / — -e p *dp = 0, * > 0, 

fl->-oo J R l 



dp 



Similar arguments are valid for the integral along T7. Thus, we have 



lim 

R— >oo 



Q (s) e s *ds 



0, t>0. 



With (j5~4")) we have that the integral over T2 becomes 



lim 

R^oo 



Q(s)e st ds 



< lim 



Q (ite 1 *) 
■ C lim / i-R tcos 



,iite 



7? 



ii?e i<#, |d^ 
= 0, t>0, 
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since cos</> < for <j> e [f? 71 "] ■ Similar arguments are valid for the integral along Tq: 



lim 

R— s-oo 



Q (s) e st ds 



0, t>0. 



Since M (s) — > c and sM (s) — >• as \s\ — > 0, ([25]) implies Q (s) « 1, as \s\ — > 0, so that the 
integration along contour L4 gives 



lim 

r->0 



Q (s) e 8t ds 



< lim 



< lim 

r->0 



Integrals along L 3 , L 5 and 7 give (t > 0) 
Q(s)e st ds = 



lim 

-R— >oo 1 1 ■ 



3 (re 1 *) 

a rtcos0 c 
1 



lim / Q(s)e st ds 

s-oo 



o 1 + (gM (ge i7r )) z 
1 



0, t> 0. 



""da, 



o 1 + (oM (ge" i7r )) 



lim / Q(s)e st ds = 2vriQ (t) . 



(55) 
(56) 
(57) 



By the same arguments as in the proof of (i) we have that (|57|) is valid if the inversion of the 
Laplace transform exists. This is true since all the singularities of Q are left from the line 7 
and appropriate estimates on Q are satisfied. Adding (|5"5")l . ([51))) and (1571) we obtain the left hand 
side of (|5"3")) and finally Q in the form given by P5"|) . 
Since 



Im 



1 



l + (gM(<ze- i7r )f 
are continuous, it follows that Q is continuous on [0, oo) . 



e~«*dq, 2Res(Q(s)e st ,s ) , i>0 



4 Example 

Suppose that F is harmonic, i.e. F (t) = Fq cos (uit) and that <p G and </> e are given by ((SJ). Then 
e and a, given by (|27|) . become 

l+as a 

£(s)=i ? 0-2- J i I 1+ 2 fc l+as° » ^( S ) = F 0^- a — 2 1+as c , S £ D. 

S Z + U Z 1 + 5 2 T , ? „ S 2 + £iT 1 + S 2 TTT^ 

In the special case a = b, which corresponds to an elastic body, we obtain 

£ 0) = F f 1 , g (s) = F 2 f 2 1 i 1 2 ' 

After inverting the Laplace transforms, we have 

Fq Fq 

s (t) — —5 cos (cot) * sin£, a (t) = — 5 — - cos (a;i) * sint, (58) 

2F . (w + l)t . (w-l)t 2F . (w + l)t . (w-l)t 
e{t) = — -sin sm , a (t) = — -sin sm . 

U! z — 1 I I LU Z — 1 2 I 
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For u> — > 1 we obtain 



£ (t) = —t sint, 



a (t) = — t sini, 



a resonance. Also, in the case when lj ~ 1 one observes the pulsation. 

We shall present several plots of e, obtained form Theorem^ (jl2"|) , in the case of the fractional 
Zener model of the viscoelastic body ©. We fix the parameters of the model: a — 0.2, b = 0.6 
and in Figures |H [5j H2 and [7J show plots of e in the cases when the forcing term is given as F = 5 
(5 is the Dirac delta distribution), F — H (H is the Heaviside function) and F (t) = cos (uit) , 
respectively. 

We see from Figure0]that the oscillations of the body are damped, since the rod is viscoelastic. 



0.5 



0.0 



-0.5 





a = 0.75 




/ a = 0.9 


| si w 






hit = 0.5^ 




\& = 0.25 




&=0.1 



Figure 4: Strain e(£) in the case F = <5 as a function of time i G (0, 38). 



The curve resembles to the curve of the damped oscillations of the linear harmonic oscillator. 
We also see that the change of a E {0.1,0.25,0.5,0.75,0.9} makes the amplitudes of the curves 
to decrease slower with the time, as a becomes smaller. This is due the the fact that a — 
1 corresponds to the standard linear viscoelastic body and a — corresponds to the elastic 
body. In the case of the forcing term given as the Heaviside function, from Figures [5] and [5] we 
observe that the body creeps to the finite value of the displacement regardless of the value of 
a e {0.1,0.25,0.5,0.75,0.9}. Creeping to the finite value of the displacement is due to the fact 
that the fractional Zener model describes the solid-like viscoelastic material. If the viscoelastic 
properties of the material are dominant (the value of a is closer to one) then the time required for 
body to reach the limiting value of strain is smaller, see Figure O compared to time in the case 
when elastic properties of the material are dominant, see Figure [S] Figure [7] shows the expected 
behavior of the body in the case of the harmonic forcing term. Namely, the oscillations of the 
body die out and the body oscillates in the phase with the harmonic function. For this plot we 
took: a = 0.45 and cj = 1.1. 

Now, we examine the case when the values of the coefficients a and b are close to each other. 
We fix them to be a = 0.58 and b = 0.6. In this case the elastic properties of the material prevail, 
since in the limiting case a = b the fractional Zener model ([JJ becomes the Hooke law for the 
arbitrary value of a G (0, 1) . We present in Figure |8] the plot of e in the case when a — 0.45 and 
oj = 1.1. In the elastic case, as it can be seen from (|58p. the frequency of the free oscillations of 
the body is u>t = 1. Since the frequency of the forcing function (a; = 1.1) is close to the frequency 
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of the free oscillations, we shall have the pulsation, as it can be observed from Figure [5] Since 
there is still some damping left (a ^ b) the amplitude of the wave-package decreases in time, as 
shown in Figure [5] 

We increase the damping effect by choosing a = 0.55, b = 0.6. The rest of the parameters 
are: a = 0.45, to — 1.1. The curve of e, presented in Figure fTOl for smaller times resembles to the 
curve of the pulsation (the elastic properties of the material prevail), while later it resembles to 
the curve of the forcing function (the viscous properties of the material prevail). 

Acknowledgement 7 This research is supported by the Serbian Ministry of Education and 
Science projects 174005 (TMA and DZ) and 174024 (SP), as well as by the Secretariat for 
Science of Vojvodina project 114 — 451 — 2167 (DZ). 

References 

[1] T. M. Atanackovic. A modified Zener model of a viscoelastic body. Continuum Mechanics 
and Thermodynamics, 14:137-148, 2002. 

[2] T. M. Atanackovic. On a distributed derivative model of a viscoelastic body. Comptes 
Rendus Mecanique, 331:687-692, 2003. 

[3] T. M. Atanackovic, M. Budincevic, and S. Pilipovic. On a fractional distributed-order 
oscillator. Journal of Physics A: Mathematical and General, 38:6703-6713, 2005. 

[4] T. M. Atanackovic, S. Konjik, Lj. Oparnica, and D. Zorica. Thermodynamical restrictions 
and wave propagation for a class of fractional order viscoelastic rods. Abstract and Applied 
Analysis, 2011:ID975694, 32 pp, 2011. 

[5] T. M. Atanackovic, S. Konjik, and S. Pilipovic. Fractionalization of constitutive equations 
in viscoelasticity. Proceedings of the Fifth Symposium on Fractional Differentiation and Its 
Applications, Hohai University, Nanjing, China, May 14-17, 2012. 

[6] T. M. Atanackovic, S. Pilipovic, and D. Zorica. Distributed-order fractional wave equation 
on a finite domain: creep and forced oscillations of a rod. Continuum Mechanics and 
Thermodynamics, 23:305-318, 2011. 



15 



e(t) 



-4 



20 



40 



eo 



100 



Figure 7: Strain e(t) in the case F(t) = cos(wi) as a function of time t £ (0, 100). 
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Figure 8: Strain e(t) in the case F(t) = cos(wi) as a function of time t £ (0, 139). 
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Figure 10: Strain s(t) in the case F(t) — cos(W) as a function of time t £ (0,600). 
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